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Abstract 

We present a scheme to accurately calculate the persistence probabilities on 
sequences of n heights above a level h from the measured n + 2 points of the height- 
height correlation function of a fluctuating interface. The calculated persistence 
probabilities compare very well with the measured persistence probabilities of a 
fluctuating phase-separated colloidal interface for the whole experimental range. 



1 Introduction 

Persistence concerns the question how long a fluctuating variable stays above a certain 
level. It is a recurrent theme in statistical physics. It started as the problem of level 
crossing in probability theory [H El |3]. If the problem is posed in the context of Gaussian 
random processes, persistence is completely determined by the correlation function of the 
stochastic variable. For these processes the essence of the problem is to calculate, from 
the correlation function, the probability that the fluctuation variable stays above a certain 
level during the time t. This turned out to be a classic unsolved problem in probability 
theory [21 HI [5] . Around the turn of the century the investigations reached a peak, with 
applications in physics, ranging from properties of the diffusion equation [6l O [8] , survival 
of spin states [9l [101 IHl IH 112] , fluctuating steps |T3] and interfaces [3 [H] , to persistence 
in order parameters [18] and non-Gaussian processes [19]. A recent review on equilibrium 
step fluctuations, which poses the problem of persistence in a wide context, is given by 
Const ant in et al. [15] . 

The investigations focussed on the calculation the asymptotic behavior. The associated 
persistence exponent sometimes characterizes a new partition of physical systems into 
universality classes. For a stationary Gaussian process the persistence probability for a 
time t decays exponentially with t. For these processes the persistence probability can 
be written as the ratio of two path integrals. The numerator involves the sum over paths 
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obeying the condition of persistence and the normahzing denominator over paths without 
the condition. For a Markov process the correlation function is an exponential and both 
path integrals can be evaluated. But as yet there is no general scheme for an arbitrary 
correlation function to calculate the persistence exponent. 

Most theoretical investigations treat the process as a continuous process in time, which 
it certainly is. However measurements of a stochastic variable are necessarily discrete in 
time. One would think that, with a sufficiently high experimental sampling rate, the 
limit of a continuous process would be seen experimentally, such that the sampling rate 
would become irrelevant. This is not obvious though, in particular for "non-smooth" 
processes A "smooth" process has a correlation function deviating from its initial 
value in a quadratic way. The influence of the discreteness of the sampling on the per- 
sistence exponent has been investigated by Majumdar et al. [16] for stationary Gaussian 
Markov processes and by Ehrhardt et al. [IT] for a number of non-Markovian smooth 
processes. In general the exponent derived from discrete sampling is lower than the con- 
tinuous exponent, because double crossing of the level in between two discrete sampling 
points are missed. These authors consider the calculation of the persistence of a discrete 
sequence more difficult than that of the underlying continuous process, as the approxima- 
tion schemes derived for continuous processes are not applicable to discrete series of data. 
For "non-smooth" processes another difficulty arises in the calculation of the persistence 
probability, because rapid fluctuations give a diverging probability on a short time scale, 
such that the mean persistence time vanishes I2U] . 

Recently we encountered the persistence problem in an experimental study of a fluctu- 
ating interface between phase-separated colloid-polymer mixtures [20]. This is an example 
of a stationary Gaussian random process, since the fluctuations are small in amplitude 
and their energy is given by a quadratic hamiltonian in the stochastic variables. On the 
time scale of the measurements the process is non-Markovian and non-smooth. We have 
collected data on the correlation function and on the joint probabilities for e.g. finding 
n successive height values above the level h. Thus here the practical persistence problem 
presents itself in a discretized form and the question arises whether one can calculate, 
from the values of the correlation function, directly the persistence probabilities of the 
same set of points in time. In an earlier publication [21] we demonstrated how to do this 
for short series. In this paper we extend the calculation to the whole set of experimental 
points. We focus here on the longer time sequences which allow for more detail than the 
similar spatial sequences. We find the values of the correlation function at equidistant 
time intervals and also the persistence probability for the same series of time steps. 

Thus experimental data were obtained in large numbers and high accuracy for an in- 
teresting system, which enable to further test the theories on the persistence problem. As 
mentioned in [T7], the continuum limit usually leads to a simplification of the calculation, 
but we find in our case the discreteness of the experimental data rather a blessing in 
disguise. Our aim is to show how the persistence probabilities for n events can be directly 
calculated from the measured n + 2 points of the correlation function for the range of 
measured points, which are neither on a short time scale nor in a fully asymptotic regime. 
We derive two sum rules for the discrete series, which play a vital role in assessing the 
accuracy of the computational scheme. 
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2 The experimental system 



Traditionally experimental studies of interfaces are carried out by means of light and 
X-ray scattering. The field obtained another dimension by experiments of Aarts et al. 
[22| 123] . in which they obtained microscopic images of fluctuating interfaces of phase 
separated colloid-polymer mixtures using confocal microscopy. Although scattering on 
interfaces is most valuable, it always yields global information on the fluctuations, while 
inspection by microscopy gives local information. However, the wave lengths and heights 
involved in the capillary waves of molecular fluids are way out of the reach of detection 
by microscopic methods. For colloidal interfaces the characteristic length and time scale 
of the fluctuations can become accessible by confocal microscopy, by lowering the surface 
tension to the nN/m range. 

Here, confocal microscopy measurements were performed on phase separated colloid- 
polymer mixtures. The colloids are 69 nm radius fluorescently labeled polymethylmetacry- 
late particles, suspended in cis/trans decalin, with polystyrene (estimated radius of gjnra- 
tion = 42 nm) added as depletant polymer. Due to a depletion induced attraction these 
mixtures phase separate at sufficiently high colloid and polymer volume fractions and a 
proper colloid to polymer aspect ratio, into a colloid-rich/polymer-poor (colloidal liquid) 
and a colloid-poor/polymer- rich (colloidal gas) phase [23]. Here the polymer concentra- 
tion acts as an inverse temperature and upon dilution the binodal is approached. 

In confocal microscopy a monochromatic laser beam is used to excite dye molecules 
attached to (in this case) the colloid. Through a dichroic mirror the outgoing light is 
separated from the incoming light. A two-dimensional confocal slice is then obtained 
through a pinhole, from single- wavelength fluorescent light emitted from the sample. For 
our experimental system, the confocal slices are only 3 colloidal diameter thick, the 
density profile between the two phases is observed as a function of fiuorescent intensity. 

A very precise location of the interface can be obtained by fitting the intensity with 
a van der Waals profile: I{z,x) = a + btanh{[z — h{x)]/c), where z is the direction 
perpendicular to the interface and x a coordinate along the interface. In the upper phase 
the density approaches a value corresponding to a + 6 and in the lower phase to a — b, 
while c measures the intrinsic width of the interface. Thus at every snapshot a function 
h{x) follows and the sequence of snapshots gives the function h{x,t). It is a practical 
separation of the particle motions, which lead at short scales to the intrinsic interface 
and the particle motions {capillary waves) which drive the long wavelengths. This opens 
up the possibility to follow in detail the motion of the height of the interface and to make 
a statistical analysis of its temporal behavior. Of course the method has its inherent 
restrictions. Just as in ordinary movie recording, snapshots have to be taken at a finite 
time intervals. For colloidal interfaces this interval can be made much smaller than the 
intrinsic time scale of the motions. 

With a Nikon E400 microscope equipped with a Nikon CI confocal scanhead, series of 
5000 snapshots of the interface were recorded at constant intervals At of 0.45 s and 0.50 
s of two statepoints to be denoted as II and IV. The latter is closer to the binodal. The 
pixels are separated by a distance Ax = 156nm and a single scan takes approximately 
0.25 s to complete. With 640 heights per snapshot, we obtain in total 640 x 5000 data 
points, which enable us to measure persistence probabilities as low as 10~^. 
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3 The correlation function 



We write the normalized correlation function as 

^7(t/g = (/i(0,t)/i(0,0))/(/i2), (1) 

where (h?) is the equal time correlation function, is the characteristic time of the 
process which we have inserted in the argument of g to make it dimensionless. It is given 
by the expression 

^ {ri + r/0 ^ ^2) 

where the ?7's are the viscosities of the two coexisting phases, Ap the density difference 
and 7 the surface tension. From the definition follows that (7(0) = 1. Heights are scaled 
with (h^Y^'^ and times with tc, such we can work exclusively with dimensionless quantities. 
Experimentally one has a discrete sampling of g{t/tc) with At as smallest interval. Thus 
one finds a sequence gn = g{ji5) with 5 = At/tc- 

Capillary wave theory for overdamped waves gives the following expression for the 
correlation function 



log(l + K^) Jo 1 + X^ 

The upper bound of the integral is k = 27r^/c? with ^ = '~f/{gAp) the correlation length 
and d the diameter of the particles. The lower bound, determined by the size of the 
system, has been set equal to 0. The upper bound is essential for the convergence of the 
integral and of influence on the short time behavior of the correlation function. Cutting 
off the capillary waves at the short-wavelength side is the poor man's way to handle the 
otherwise diverging interface width {h'^). Fig. [1] shows the function (|3]) together with data 
points referring to two experiments, denoted by II and IV [20j. The choice of 6 is actually 
a fit of tc- So we have plotted in the figure the data points with t/tc = n6 and 6 = 0.04 
for II and 6 = 0.02 for IV. 

Although the data points follow the curve reasonably well, the scatter of the data 
points is manifest and more important, though only visible in the insert, is the fact that 
the first few points are well below the curve. In this paper it is of less importance how 
well the experimental correlation function can be represented by a theoretical curve, since 
we are interested in the problem to directly calculate the persistence probabilities from 
the experimental correlation function. 



4 Persistence probabilities and sum rules 

To formulate the persistence probability for a discrete series, we form the n x n matrix 
gij = g\i-j\, which is determined by the first n values of gn- Using that the process is 
Gaussian, the probability on a sequence of values {hi, ■ ■ ■ , hn) is given by the formula 

G{hi, ■■■,hn) = exp f-i Ji,jhih}\ , (4) 
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Figure 1: Height-height correlation functions, curve ([3]) and experimental points; circles 
refer to 6 = 0.04 (II) and squares to 6 = 0.02 (IV). Insert: first 10 measured correlation 
points and curve ([3]). 

with the normalization D = (27?)" det where det g is the determinant of g. Jij is the 
inverse of gij, which is again a. n x n matrix and its matrix elements also individually 
depend on n (while those of g only depend on \i — The justification of @ is based on 
the fact that the correlator (hihj) as calculated with @ leads indeed to g\i-j\- 

The expression for the persistence probability p„(/i), on a sequence of precisely n events 
above h, follows as the ratio of two integrals 

Pn{h) = g-("+)-(/i)/g-+(/i). (5) 

The numerator is the (n + 2)-fold integral 

q-(n+)-^h)= I dh,G{ho,---,hn+i), (6) 

where the integration domain P, indicated by the superscript on q, is given by the condi- 
tions h > Hq > — oo, h < hi < CO, and h > hn+i > — oo. It selects the sequence of events: 
starting with a value below h, followed by n points above h and terminated by a value 
below h, which are the sequences of precisely n successive values of the height hi > h. 
Similarly the denominator is given by the double integral 

/h roo 
dho dh^G{ho,hi) (7) 
-oo J h 
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The inverse matrix Jj j depends for n = 2 only on the first value gi of the correlation 
function. So the q ^{h) is generally given by 

n-+(h) - r Hh r Hh ^M-{hl + hl-gihoh)/{l-gj)] 
q {h)-J_JhoJ^ dh [2n{l-g!)]y^ 

The denominator serves as normalization, since one has the relation 



g-+(/,) = ^g"(«+)-(/,) (9) 
n=l 



which gives the sum rule for the total probability 

EPnW = l. (10) 

n=l 

The derivation of is based on repeated expansion of the relation 

q-+(h) = q-+-{h) + q-++{h) = q'^' {h) + q-++-{h) + = ■ ■ ■ , (11) 

which simply states that the probability to find a sequence — h is the same as finding it 
from the events — h + and — I — . Together they extend the integration over the last 
height variable over all values and the probability on — h results. Repeatedly adding a 
new point in the sequence ending with + leads to the identity (jH]). 

There exist one other sum rule. Consider q'^{h), which is the single integral over Hq 
with ho > h. The expansion, similar to ( ITTj) starts as 

q+{h) = q+~{h) + q++{h) = q-+-{h) + q++- + g"++(/i) + g+++(/i). (12) 

Systematically replacing every + at the begin or end of the string by the sum of strings 
extended with a + and — gives 



g+(/i) = ^ ng-("+)-(/i). (13) 

n=l 



This yields for the mean value the sum rule 

Y.^Pn{h) = q+{h)/q-+{h). (14) 

n=l 

We have not found further sum rules. 

The sum rules f llOp and fll4p . which generally apply for stationary Gaussian processes, 
involve doable integrals. q~^{h) is an error function and q ^{h) can be reduced from the 
double integral ([8]) to a single integral, for which an analytic expression exist for h = 

m 

g-+(0) = - - - arctan | ) . (15) 

The sum rules are typical for the discrete series as the continuum limit for non-smooth 
processes is singular. If one could straightforwardly define in the continuum limit a 
probability density p{t) on a persistence interval t, one would expect that it would be 
related to the discrete Pn as 

p{nAt)At ~ Pn, (16) 
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where we suppressed the h dependence for the moment. The probabihty density is then 
normahzed ^ 

/ rftp(t)~^p„ = l. (17) 

Consequently one would expect the mean of the discrete series to diverge as 

^np„ = (At)-i^(nAt)p„~(At)-W dttp{t) ^ {At)-\ (18) 

n=l n=l ''^ 

However the expression (flSjl shows, with ~ 1 — 0{At), that 

g-+ ~ (A^)V^ (19) 

while g"*" is independent of At. So the discrete mean diverges as (At)~^/^ in contrast to 
the expected behavior f fTSj) . The discrepancy is due to the Brownian fluctuations on all 
length scales for a non-smooth correlation function^ |2Q]. This gives a non-integrable 
probability density for short times, making the normalization f lT7|) questionable. 

We calculate the Pn{h) and then check whether total and mean corresponds with the 
exactly calculable ratios. In particular f[T^ is a stringent test on the the calculated values 
of p„(/i), since it more sensitive to the large n values than (fTOj) . 

The restricted integration domain prevents the integral (j6]) from straightforward eval- 
uation. Although it involves a finite set of integrations, which indeed can be performed 
by standard techniques for small n, a direct evaluation of (E]) is impossible for large n. As 
we shall show, expression (jS]) has a definite asymptotic large n behavior, but the exper- 
imental data are not at all exclusively determined by this asymptotic behavior. In fact, 
the practical regime for which accurate data can be collected, shows important transient 
behavior. 



5 The Markovian case 

Our calculational scheme is inspired by the perturbation technique of Majumdar and Sire 
[1], which takes the Markovian case as lowest approximation. The Markovian case has an 
exponential decaying correlator 

(?(t/g =exp(-At/g, (20) 

such that Qn = (?" with g = exp(— A5). The inverse matrix J is then a band matrix with 
all elements except on the diagonal: 

JqO = Jn+l,n+l = 1/(1 ^ fi'^)' Jii = (1 + 9^)^00^ Ji,i±l = ~gJoO- (21) 

Note that this also holds for a matrix of finite dimension n. If the matrix J is restricted to 
the diagonal and the subdiagonals, the corresponding joint probability can be factorized 
in several ways. We present here the symmetric representation for the above Markovian 
case ^ 

M(/io, ■ ■ ■ , K+i) = fo{ho) (U K{h,_u hi)^ Mhn+i). (22) 

The initial (and final) function is given by 

/o(x) = exp[-«oX^]/[2vr]V2 (23) 



7 



with Mo = 1/4. The kernel reads 



K(x,y) = '^^-f^^^:^ , (24) 



with the values 

" = ^(137)' -=(Y37)- (25) 

Physically, the factorization results from the fact that in a Markovian process the proba- 
bility on an event only depends on the probability of the previous event. Mathematically, 
any matrix which is restricted to the diagonal and the subdiagonals, can be seen as a 
Markovian matrix. We use the more convenient symmetric form, allowed by time re- 
versal symmetry, rather than the standard conditional probability, with an asymmetric 
kernel. We will lean heavily on choosing the optimal Markovian approximation, which 
uses optimal values for the u and v and not those tied in with g. 

Any of these representations give the values of Pn{h) recursively. Define a set of 
functions /„(a;) with fo{x) given by (1231) . /i is constructed from /o as 

h{y)= r dxfo{x)K{x,y). (26) 

J —oo 

The general term is defined by recursion for 1 < i < n + 1 

POO 

My)= / dxf,.i{x)K{x,y). (27) 

J h 

Then Pn{h) can be expressed as 

pM = dxU^,{x)fo{x)/q~+{h). (28) 

J —oo 

In the Markovian case the multiple integral (Q becomes a repeated integral transformation 
with K as kernel. Asymptotically the result is dominated by the largest eigenvalue of the 
kernel, which then leads to a persistence exponent. 



6 Calculational scheme 

Our approximation scheme is based on a separation of J in a Markovian part M and a 
remainder H 

J = M+n. (29) 

is a Markovian approximant, i.e. a matrix which is confined to the diagonal and the 
subdiagonals. Thus one can relate a kernel K to A4 as in ( |22|) . but with as yet unspecified 
values of Uo,u and v. Also a Markovian joint probability M can be associated with this 
part (as in (|22l) . which enables to define a Markovian approximation Pn{h) 

Pl{h) = [ dh^M{ho, ■ ■ ■ , hn+i), (30) 

q ^{n) Jv 

with V the same integration domain as in ([6]). The full probability Pn{h) then reads 

Pnih)=pl{h){exp{-n))o. (31) 
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where the average {A)o is defined as 



/-P dhj M{ho, ■■■ , hn+i)A{ho, ■ ■ ■ , 
/p dhi M{ho, ■ ■ ■ , hn+i) 




(32) 



The exponential in fl3T|) will be evaluated by the cumulant expansion. The Markovian 
average of the first cumulant of "H as well as the higher cumulants are calculated in a 
similar iterative way as Pn{h). 

There is a large freedom in choosing the Markovian part. Any matrix with only a 
non-zero diagonal and subdiagonals would do. This freedom can be exploited through 
the inequality 



The optimal Markovian matrix gives the largest value for the right-hand side. The cal- 
culation thus consists of finding an optimal Markovian approximation and evaluating the 
corrections. For the convergence of the calculational scheme the optimization is crucial. 
As optimization parameters we can use Uq and for each kernel a choice for u and v. The 
practical optimization is a trade-off between the optimum and calculational simplicity. It 
is expedient to have all kernels the same, with the exception of the first and last. They 
have anyway a different role, carrying the system from below h to above h, while the 
other progagators keep the system above h. So we remain with an initial (and final) Uq, 
the first (and last) pair ui, Vi and the bulk pair u, v. A further restriction stems from the 
chosen integration procedure. The integrals can be fast and accurately evaluated using 
Gaussian quadrature, with a Gaussian as weight factor. The initial /o(a;) is a Gaussian. 
Requiring that after an iteration the asymptotic behavior is still the same Gaussian, u, v 
and Uq have to be related as = f ^ + Mq. A constant Gaussian allows to use the same 
positions and weights in the Gaussian quadrature for each iteration. This reduces the 
freedom from 5 to 3 parameters: Uq,Vi and v. 

Thus our calculational scheme starts out by finding for every n the optimal values 
of the parameters uq, vi and v. According to the inequality (|33|) this gives persistence 
probabilities Pn{h) which are too low. Consequently the total persistence and the mean 
as calculated with this approximation will be lower than the exact values calculated from 
the sum rules ffTU]) and f lT^ . Using the the Gaussian quadrature this is a fast routine. 
The next step is the calculation of the second cumulant, which leads to an overestimation 
of the persistence probability. Thus including the second cumulant the sum rules are 
approached from above. As the second cumulant requires to calculate many correlation 
functions (going up with n^) this becomes already a longer calculation. It takes about an 
hour on a simple PC to calculate persistence probabilities up to n = 150. Needless to say 
that the computation of the third cumulant, where the number of correlation functions 
goes up as n^, is even more time consuming. In practice we could calculate the third 
cumulant up to n = 40 in a reasonable amount of time. Fortunately the changes in the 
persistence probabilities due to the inclusion of the third cumulant are so small that they 
do not infiuence the value of the sum rules. 

7 Test of the approximation scheme 

Before we carry out the computation of the persistence probabilities of the actual data, we 
inspect a clean case by taking the correlations points gn from the analytic expression (EJ. 



{exp{-n))>expi-{n)). 



(33) 
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This has the advantage of noise- free input and allows to vary the sampling rate 5. The 
test consists out of the calculation of the sum rules with the cumulant expansion outlined 
in the previous section. We consider two samplings of the curve ([3]), 5 = 0.04 and twice 
as narrow 6 = 0.02, which correspond closely to the experiments II viz. IV [21]). We 
take K = 50 in both cases. In Table 1 we have summarized the sum rules for the total 



sum 


II: 5 = 0.04 


IV: 5 = 0.02 


1st 


2nd 


exact 


1st 


2nd 


exact 


EnPn(O) 


0.9935 
4.6998 


1.00002 
4.8517 


1.0000 
4.8420 


0.9932 
6.3623 


1.0005 
6.5417 


1.0000 
6.5351 



Table 1: Sum rules for the testcurve; cumulants and exact value. 



f lTU]) and and mean f ll4l) persistence in the 1st and 2nd cumulant approximation, for the 
two samplings of ([3]) with level h = as discriminator. We have calculated 100 points for 
sampling II and 150 points for IV, which suffice to saturate the values of the sums. The 
third entry, exact, gives the values of the sum rules calculated from the ratios of integrals. 
The accuracy in the total probability according to the first cumulant, giving a rigorous 
lower bound, is amazing. For the mean value this is still impressive, but indicates that 
the probability for the higher values of n is somewhat too low. The second cumulant, 
always overestimating the probabilities, makes well up for this deficiency of the lower 
bound. The overall impression is that the scheme performs very well for points taken 
from the curve We have doubled the cut-off k to see the infiuence. It marginally 
lowers the curve, except that gi is decreasing visibly. A larger cut-off makes the process 
less "smooth". This gives a deterioration of the bound. We have verified that including 
the third cumulant does not change the numbers in Table 1 appreciably. 

8 Comparison with the experiments 

The challenge is whether the calculation of the persistence probabilities from the data 
points of the correlation function is also accurate. First we test this again with inspecting 
the sum rules, which are summarized in Table 2. The entry "exact" refers to the 
calculation of the sums using f llOp and fll4p . which requires only the value of gi. The 
marks around "exact" refer to the uncertainties in gi determining the ratios, which is 
not correlated to the uncertainty in the other gn, together determining the persistence 
probabilities. The sums for statepoint II are based on 100 points, which is about the 
number of reliable values of gn- The sums of state point IV are extended to 150 points. 
In both cases the value of Pn(0) is then so small that inclusion of further points will not 
change the sums. The lower bound is considerably less accurate than for the generated 
gn- This is not surprising. The above mentioned scatter in the data adds to the non- 
Markovian character. Even if the process were strictly Markovian, the noise would be 
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sum 


state point II 


state point IV 


1st 


2nd 


"exact" 


1st 


2nd 


"exact" 


E„np„(0) 


0.9480 
2.8760 


1.0024 
3.7093 


1.0000 
3.6513 


0.9339 
3.2761 


1.0034 
4.9344 


1.0000 
4.7933 



Table 2: Sum rules for the experiments; cumulants and calculated mean. 



seen as a deviation from Markovian behavior and would lead to a convergence question. 
The second reason is that the measured gi is quite low: gi = 0.65 (II) and gi = 0.79 (IV). 
We noticed already that less "smooth" processes have a slower convergence. However the 
second cumulant approaches nicely the correct values, indicating that these calculated 
persistence probabilities are accurate. 

The real challenge is to see how well the persistence probabilities calculated from the 
measured correlation points compare with measured probabilities. This is shown Fig. [2] 
for the two state points (with h = 0). The calculated points are taken from the second 
cumulant approximation, since the third cumulant has no appreciable influence. The 
agreement for state point IV is as good as one can hope for probabilities as small as 10~^. 
For state point II there is a systematic deviation for larger n. This may be caused by 
artifacts resulting from the confocal slicing, (/i^) and therefore the resolution relative to 
{h"^) are higher for IV than for II, leading to more pronounced artifacts in II [20]. Note 
that the asymptotic exponential decay sets in around n = 30 which takes as long as 15s 
in real time. 

9 Discussion 

We have presented measurements for the correlation function and the persistence proba- 
bility of the fluctuating heights of a colloidal interface. We have developed a calculational 
scheme which enables to find the persistence probability directly from the measured cor- 
relation function. This means that our scheme is independent of the agreement of the 
correlation function to a function like ([3]). The only assumption is that the distribu- 
tion (jl]) is Gaussian, which stems from the fact that the thermal interface fluctuations 
are small deviations from equilibrium. The predictions of the persistence probabilities 
from the measured correlation function agree very well with the measured persistence 
probabilities. 

Finally we like to make a few comments on the results. 

1. Fig. [2] shows that the curve approaches an exponential decay, albeit that a window 
from 30-100 points is not large for an accurate determination of the exponent. As we 
mentioned the two experiments may be seen as two samplings of the same correlation 
function. In order to compare the two samplings we translate the behavior to the 
time domain. For the two samplings we estimate the decay of Pn(0) ~ exp{—ilj{6)n). 
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Figure 2: Probabilities (lines) calculated from the measured correlation function and 
measured probabilities (points) for state point II (lower curve, squares) and IV (upper 
curve, circles) and h = 0. 



Since n6 = t/tc, the persistence exponent 9 equals 9 = ip{5)/5. For the smallest 
5 we find 9 = 0.027 and for the larger 6,9 = 0.026. These two values need not 
be the same. As observed by Ehrhardt et al. [17], the exponent for non-smooth 
processes is quite sensitive to the sampling rate. The discrete series does not record 
what happens in between the measured points. The continuous process might dive 
below the level h between recorded points. Such paths are excluded in a continuous 
formulation, but are included in the discrete version [161 120] • This leads to a smaller 
exponent for a larger 6. 

2. We note that we cannot take the continuum limit for practical and essential reasons. 
Making our 6 smaller requires longer series to calculate as the two samplings show. 
But more importantly, a smaller 6 will lead to a sharper peaked kernel, which 
ultimately has to approach a 5-function. Such a delta peak cannot be treated by 
Gaussian quadrature, which has helped to speed up the integrations by about a 
factor thousand with respect to e.g. a Simpson rule. 

3. We have checked that higher cumulants have virtually no influence on the persistence 
probabilities calculated here. We noticed, however, that although the third cumulant 
is small, it develops a linear dependence on n. It has no influence on the presented 
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sum- rule results nor on the presented persistence probabilities, since it changes only 
the probabilities which are already too small to contribute. But it will have a non- 
negligible influence on the decay exponent of the discrete series, showing up in a 
region which is beyond our measurements. This also means that the experiments 
have not entered the fully developed asymptotic regime. 

Acknowledgement The authors thank Wim van Saarloos and Henk Hilhorst for 
stimulating discussions. 



References 

[1] D. Slepian, Bell Syst. Tech. 41 (1962) 463. 

[2] M. Kac, SIAM (Soc. Ind. Appl. Math.) Rev. 4 (1962) 1. 

[3] I. F. Blake and W. C. Lindsey IEEE Trans. Info. Theory 19 (1973) 295. 

[4] S. N. Majumdar and C. Sire, Phys. Rev. Lett. 77 (1996) 1420. 

[5] J. Krug, H. Kallabis, S. N. Majumdar, S. J. Cornell, A. J. Bray and C. Sire, 
Phys. Rev. E 56 (1997) 2702. 

[6] S. N. Majumdar, C. Sire, A. J. Bray and S. J. CorneU, Phys. Rev. Lett. 77 (1996) 
2867. 

[7] T. J. Newman and Z. Torockzai, Phys. Rev. E 58 (1998) R2685. 
[8] H. J. Hilhorst, Physica A 277 (2000) 124. 

[9] B. Derrida, A. J. Bray and C. Godreche, J. Phys. A 27 (1994) L37. 

[10] B. Derrida, V. Hakim and V. Pasquier, Phys. Rev. Lett. 75 (1995) 751. 

[11] B. Derrida, V. Hakim and R. Zeitak, Phys. Rev. Lett. 77 (1996) 2871. 

[12] S. N. Majumdar and A. J. Bray, Phys. Rev. Lett. 81 (1998) 2626. 

[13] D. B. Dougherty, I. Lyubinetsky, E. D. Wiliams, M. Constantin, C. Dasgupta and 
S. Das Sarma, Phys. Rev. Lett. 89 (2002) 136102. 

[14] H. Kallabis and J. Krug, Europhysics Lett. 45 (1999) 20. 

[15] M. Constantin, C. Dasgupta, S. Das Sarma, D. B. Dougherty and E. D. Williams, J. 
Stat. Mech. (2007) P07011. 

[16] S. N. Majumdar, A. J. Bray and G. C. M. A Ehrhardt, Phys. Rev. E, 64 (2001) 
015101. 

[17] G. C. M. A Ehrhardt, A. J. Bray and S. N. Majumdar, Phys. Rev. E, 65 (2002) 
041102. 



13 



[18] K. Oerding, S. J. Cornell and A. J. Bray, Phys. Rev. E 56 (1997) R25. 

[19] O. Deloubriere and H. J. Hilhorst, J. Phys. A. Math. Gen. 33 (2000) 1993. 

[20] V. W. A. de Villeneuve, J. M. J. van Lccuwcn, J. W. J. de Folter, D. G. A. L. Aarts, 
W. van Saarloos and H. N. W. Lekkerkerker EPL 81 (2008) 60004. 

[21] V. W. A. de Villeneuve, J. M. J. van Leeuwen, W. van Saarloos and H. N. W. Lekkerk- 
erker J. Chem. Phys. 129 (2008) 164710. 

[22] D. G. A. L. Aarts, M. Schmidt and H. N. W. Lekkerkerker Science 304 (2004) 847. 

[23] C. P. Royal, D. G. A. L. Aarts and H. Tanaka, Nature Physics, 3 (9) (2007) 636. 

[24] H. N. W. Lekkerkerker, W. C. K. Poon, P. N. Pusey, A. Stroobants and P. B. Warren, 
Phase behavior of colloid + polymer mixtures. Europhys. Lett. 20 (1992) 559. 

[25] U. Jeng, S. Esibov, L. Crow and A. Steyerl, J. Phys.: Condens. Matter, 10 (1998) 
4955. 



14 



